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TWO SUBROUTINES FOR 8K FORTRAN 


DECUS Program Library Write-up 


DECUS NO. 8-480a 


INPUT 

A reloctable input routing for input in free format to 8K fortran programs 

Usage: 

CALL INPUT (IDEV, A, NR) 
where 

IDEV is device code for input (standard 8K fortran usage (device 
codes 1 and 2 only), for device code 0 see below) 

A is first position to store the incoming numbers in (Floating pt 
vector) 

NR is number of items to be input 

Input is in a free format, the numbers are separated (and terminated) by , 
or return. Two terminators in sequence without intervening input will set 
number®0. 

Special operators are as follows: 


Character 

ASCII 

Meaning 

9 

254 

terminates number 

cr 

215 

" " must be given before return to 

main program is effected (but 
see CTRL/Z) 

e 

256 

decimal point 

rub 

577 

removes last input before terminator is given 

ctrl/u 

f 

225 

556 

removes input since last cr 

restarts input routine (removes all input) 

ctrl/z * 

241 

terminates input & returns to main program, all 
remaining numbers are set to 0 

- 

255 

negative number 

ctrl/c * 

205 

return monitor directly. 


All other characters and illegal use3 of the above (or excess figures after 
all numbers have been input) are ignored and not echoed on the 
teletype. 


#t Z, and -f C will function during high speed reader input. 

Ex CALL INPUT (i, A(5), 10) will read from the teletype 10 numbers to be 
stored in A(5)-- A(l5;. 

A call to input with device number 0 will set the referenced variables = 0. 
This is a fast and echonomical (in core space) method of zeroing a vector 
before use 

e g CALL INPUT (0,A(5), 10) 

will set A(5) - A(15) to zero using less core space than the equivalent do 
loop. 
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The routine requires 4 pages but much core space is saved by omitting read 
statements and their FORMATS not to mention the simplicity of input as 
opposed to formated input. 

The routine ignores the parity bit and can be used on all teletypes. If the 
routine is to be used in a non-PSO environment bo must some slight re- 
editing be done, change CALL 0,EXIT at location »fPEXIT" to "HLT" or "RETURN 
INPUT" whichever suits you best. 


RANDU and GAUSS, random number generators for 8K fortran, 


Usage. 

RANDUi CALL RANDU (A, X) 

GAUSSi CALL GAUSS (AM, S, A, XN) 

Parameters i 

~ 4 n 0n f f rSt °f llin S be contain any random floating point 

fn r r t iT n U c ? ntains a new f *P‘ number which should 

sho^d Xt ahould not Or Caution A 

should not be manipulated by the SABR rutines for f.p. numbers 

the C nmb£r! ° rifflnal definition ) as these routines normalize 

On return a random normal number between 0 and 1 
Mean of required population 
Standard deviation of required population 
random hormal number from the population (M, S) 


X 

AM 

S 

XN 


Method* 

DSC “ S 5 ‘ 25 written in SABR. Th. final maker ia 
roJ thi Jairto. “ , P K Si e floating point number. The 36bit starting number 

(".! nnlT.T^nSSS)?" “ 3t ° re,i ^ * “ d * ranSferred ’ , ° r4 b > 

SJthode^or^oi^jnf^a 0 " *? n ° r "““ 1 “otriOntions (from RW Hamming Numerical 
Methods for scientists and engineers, McGrawhill 1962). * J ’ 
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K 

Z X. - K 
i - 1 1 2 

Y " — ’ - 

fk -12 

with K « 12 this reduces to 

12 

Y * Z X. -6 

i * 1 1 

and adjusted to the required M and SD we have 

XN = Y*S + AM 

The are generated by 12 calls to RANDU 

Execution time (total, including transfer from calling fortran program) 

RANDU approxt 30 msec 
GAUSS •' 150 msec 

Storage 3 pages 
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